Geometric stochastic resonance in a double cavity 
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Geometric stochastic resonance of particles diffusing across a porous membrane subject to oscil- 
lating forces is characterized as a synchronization process. Noninteracting particle currents through 
a symmetric membrane pore are driven either perpendicular or parallel to the membrane, whereas 
harmonic-mixing spectral current components are generated by the combined action of perpendic- 
ular and parallel drives. In view of potential applications to the transport of colloids and biological 
molecules through narrow pores, we also consider the role of particle repulsion as a controlling factor. 

PACS numbers: 05.40.-a, 05.10.Gg 



I. INTRODUCTION 

Historically, research on Stochastic Resonance (SR) p] 
focused mostly on systems with purely energetic poten- 
tials, either continuous or discrete. However, as pointed 
out in Ref . [2] , in soft condensed matter and in a variety 
of biological systems [21 0], particles are often confined 
to constrained geometries, such as interstices, pores, or 
channels, whose size and shape can affect the SR mech- 
anism [5 . Indeed, smooth confining geometries can be 
modeled as entropic (i.e., noise or temperature depen- 
dent) potentials [5], capable of influencing the response 
of the system to an external driving force (see, for a re- 
view, Ref. [7]). 

A more interesting scenario emerges in the case of 
sharp confining geometries. In a recent paper [5] we 
showed that a Brownian particle confined to two distinct 
cavities divided by a porous medium, say a membrane, 
undergos SR when driven by an ac force perpendicular to 
the membrane. This means that, at variance with ordi- 
nary SR PP, optimal synchronization between drive and 
particle oscillations for an appropriate noise level, occurs 
even in the absence of a bistable effective potential, of 
either energetic or entropic nature However, such a 
manifestation of SR in higher dimensions requires adopt- 
ing extremely sharp geometrical constrictions to separate 
the two cavities. The magnitude and conditions of such 
a geometric SR effect are sensitive to the geometry of the 
cavities and the structure of the pores. 

This instance geometric SR should not be mistaken for 
the socalled entropic SR introduced in Ref. [5j , where the 
Brownian particle can switch cavity only by overcoming 
the entropic barrier determined by the geometric con- 
striction associated with a smooth pore. In the absence 
of an energetic barrier, the entropic barrier alone deter- 
mines the magnitude of the SR effect that occurs when 
a periodic force drives the particle across the pore, al- 
though the evidence reported by Burada et al. [5] hints 
at an interplay of entropic and energetic barriers, rather 
than to a mere entropic effect. 

Geometric SR is a peculiar manifestation of driven 



Brownian motion in septate channels [9HT3|. In these 
channels compartments are separated by zero-thickness 
partition walls and the pores are modeled by structure- 
less holes pierced at the center of the partition walls. A 
septate channel cannot be analyzed in terms of the re- 
ductionist approach of Ref. [5], i.e., as entropic channels 
[H [7] , because the geometry of its pores turns out to be 
too sharp for a one dimensional (ID) kinetic equation 
approximation to hold |14) . However, sharp pore geome- 
tries are known to enhance most of the noise controlled 
transport mechanisms proposed in recent years pQ, thus 
making this class of channels particularly suitable for ex- 
perimental demonstrations. 

In this paper we restrict our analysis to two dimen- 
sional (2D) geometries. Our choice is partly motivated by 
the growing interest in vortex superconducting devices, 
a class of artificial devices with numerous technological 
applications, including flux qubits, SQUIDs and super- 
conducting rf filters [121 [H] ■ Superconducting samples 
with two vortex boxes connected by a thin pore of al- 
most any geometry can be easily fabricated. Vortices are 
trapped inside the boxes with binding energy of the order 
of QgLt/X 2 , where <l>o is the magnetic flux quantum, A 
is the London penetration depth, and L t is the depth of 
the two vortex traps. Magnetic vortices repel one another 
through a logarithmic pair potential, while their density, 
n = H/&o, is controlled by the intensity H of the ap- 
plied magnetic field. In the dilute limit, H < $o/A 2 , the 
vortex-vortex interactions become negligible, so that the 
transport properties of a single trapped vortex are not 
overshadowed by many-body effects, ac drives and noise 
sources can be easily implemented as Lorentz forces gen- 
erated by independent electric currents injected into the 
sample parallel and perpendicular to the pore axis. De- 
tection of SR under such experimental conditions is thus 
regulated by the applied current sources only; in partic- 
ular, the noise parameter can be varied independently of 
the constant operating sample temperature. 

This paper is organized as follows. In Sec. [ll]we sum- 
marize the spectral properties of geometric SR, as first 
reported in Ref. [8]. In Sec. Ill we analyze geomet- 



2 




-q 1 — i — i i 1 1 1 1 1| 1 — i — i i i 1 1 ■ ■ | 1 — i — i i i 1 1 1 ■ | 1 — i — i i i 1 1 1 ■ | 

10" 5 10" 3 D 10" 1 



FIG. 1: (Color online) (a) Brownian particle confined to a 
2D box divided in two compartments by a partition with an 
opening of width A at the center, (b) Geometric stochas- 
tic resonance. x(D) versus D for different values of Q at 
Ao — 0.045 (main panel) and A x at Q. — 0.01 (inset). Other 
parameters are: xl = ul = 1 and A = 0.1. In the inset 
we also report our predictions for the SR peak position D max 
(vertical arrows) and height x(D max ) (top line). In both plots 
the dashed line represents the decay law x(D — > oo) (see Sec. 

m. 



ric SR as a noise controlled synchronization mechanism. 
The distribution of the crossing times of a single particle 
through a pore is investigated in the case of both longi- 
tudinal and transverse ac drives. In Sec. |IV| we study 
the particle current through a pore under the combined 
action of longitudinal and transverse drives with incom- 
mensurate frequencies. Harmonic mixing is detected as a 
consequence of confinement, even in the absence of non- 
linear particle interactions. In Sec. [V]we explore the role 



of particle repulsion in controlling geometric SR, by vary- 
ing the particle density and the pair interaction intensity. 
Finally, in Sec. |VI| we add some concluding remarks. 



II. GEOMETRIC STOCHASTIC RESONANCE 

Let us consider an overdamoed Brownian particle 
freely diffusing in a 2D suspension fluid contained in two 
symmetric xl X Ul compartments, with reflecting walls 
|17) . connected by a narrow pore of width A, as illus- 
trated in Fig. [ija). The overdamped dynamics of the 
particle is modeled by the Langevin equation, 



dr 
dt 



= -A(t) +VD£(t), 



(1) 



where A = (A x ,A y ) are the x, y components of the 

driving force and £(i) = {£, x (t), £, y (t)) are zero-mean, 
white Gaussian noises with autocorrelation functions 
(&(*)&(*')) = 2<%£(f - f), with i,j = x, y. Equation {l} 
has been numerically integrated by a Milstein algorithm 
[18] . Stochastic averages were obtained as ensemble aver- 
ages over 10 6 trajectories with random initial conditions; 
transient effects were estimated and subtracted. 

In the presence of a longitudinal ac drive oriented 
along the x axis, A x {t) — A x cos(fii), the Brownian 
trajectories embed a persistent harmonic component, 
x(D) cos[ili — <fr(D)], whose amplitude, x is plotted versus 
D in Fig. |TJb) . For details about the curves 4>(D), not 
relevant to the discussion below, the reader is referred to 
the original paper [5]. 

The occurrence of geometric SR is clearly documented 
in Ref. [8j- Its main distinctive properties that can be 
summarized as follows: 

(i) x(D) peaks for an appropriate noise intensity, D mstx , 
with upper bound, 



x(D 

max J _ ^L- 
7T 



(2) 



A satisfactory estimate of -D max was obtained by match- 
ing the half-drive period with mean first-exit time from 
one compartment at A x = 0, ta(D), that is, 



T A(-Dmax) 



To 



7T 



(3) 



Note that ta(-D) is inversely proportional to D and 
strongly depends on the pore width, A, as computed in 
Rcf. [12] . Hence, D max is proportional to fi; 

(ii) SR is restricted to Ao > A c or f2 < f2 c , as a conse- 
quence of the geometric condition, 



A x A 

il 7T 



(4) 



This is an important difference with respect to ordinary 
SR pQ, where there exist no such onset thresholds in the 
drive parameters space; 
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(iii) x(D) obeys the approximate SR curve, 

x (D) 



where, 



x(D) 



x {D) 



V/I + ^TA^)] 2 ' 

A x x L /D _ 
tanh (A x xl/D) 



(5) 



(6) 



is the amplitude of the (x(t)) oscillations in the adiabatic 
limit f2 — > 0. As a consequence, for large noise, x{D) oc 
1/D. 
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FIG. 2: (Color online) Synchronization mechanism. Distri- 
bution density of the residence times for different periods 
Tq — 2n/Q (in the legends) of the (a) longitudinal drive with 
A x = 0.05 and (b) transverse drive with A y — 0.5. The cav- 
ity dimensions are as in Fig. [TJand D — 0.015. Inset in (a): 
height of the first peak, n — 1, versus Q. Vertical arrow de- 
notes the resonant f2 value predicted in Eq. Q with ta given 
in Q2]. 



The properties listed above follow from simple geo- 
metrical considerations. Equation ([3| is the geometric 
counterpart of the standard SR condition [T], where, in 
the absence of an energetic barrier, the Arrhcnius time 
is replaced by an appropriate diffusion time across the 
pore. As is apparent in Eq. (JTjj), for weak noise and low 
drive frequencies, A x (t) presses the particle against the 
walls of the container opposite to the dividing wall and, 
as a result, the average particle displacement, (x(t)), ap- 
proaches a square waveform with amplitude x^. The 
amplitude of the Fourier component of (x(t)) with an- 
gular frequency f2 is then, Axl/k. From this remark 
it immediately follow the inequalities in Eqs. ^ and 
Q, A x /fl being the dr iven oscillation amplitude of an 
unconstrained Brownian particle. The upper bound of 
Eq. ([2]) hold for vanishingly low f2, see Fig. [ljb). For 
stronger noise but still low drive frequencies, xo[D) tends 
to A x t. x , where t x = x 2 L /iD is the longitudinal diffusion 
time across a cavity compartment. The pore effectively 
suppresses the particle oscillations, with damping con- 
stant ta(-D) -1 , only at relatively high drive frequencies, 
as apparent from Eq. ([5]). The consistency of these an- 
alytical results with the simulations is quite satisfactory, 
as shown in [5]. 



III. SYNCHRONIZATION MECHANISMS 

As summarized in Sec. [TTJ the evidence for geometric 
SR reported in Ref. [5] focuses on the D dependence of 
the harmonic component of (x(t)} with driving frequency 
fl. This spectral characterization of geometric SR will be 
further discussed in the following section. 

An alternative and perhaps deeper insight into the un- 
derlying resonance mechanism can be gained by consider- 
ing the so-called synchronization characterization of SR 
[TU] . In this approach one looks at the residence times, T, 
of the Brownian particle in either cavity - which one is ir- 
relevant, due to the mirror symmetry of the process with 
respect to the cavity divide. Their distribution densities, 
7V(T), exhibit a prominent peak structure and, more re- 
markably, a resonating D and T dependence. 



A. Longitudinal drive 



Let us consider first the case of longitudinal drives, 
A„ = 0. 



In Fig. [§ £ 



we plotted N(T) for different values 
of fi at constant D. It is immediately apparent that the 
N(T) peaks are centered around 



2 1 T " 



1,2,3, 



(7) 



due to the fact that the left-to-right pore crossings are 
most likely to occur when the longitudinal ac force A x (t) 
points to the right and vice versa. Such a synchroniza- 
tion mechanism has been discussed at length in the SR 
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literature pQ. A detailed numerical analysis (shown in 
Fig. 2(a) for n = 1) proves that the height of (or, more 
precisely, the area enclosed under) the n-th distribution 
peak increases with CI up to an optimal value Cl = Cl n , 
and then decreases for higher Cl. Note that Cl n is deter- 
mined by the optimal synchronization condition |19j 



n--j T n =T A (D), n= 1,2,3,.... (8) 

Note that for n = 1 this equality coincides with the spec- 
tral SR condition in Eq. pj). This proves that the first 
peak (but this conclusion actually applies to all them!) 
can be enhanced to a maximum by acting either upon Tq 
or ta- Hence, geometric SR is also a bona-fide resonance 

US]. 



B. Transverse drive 

We now consider the case of transverse drives, A x = 0. 
Clearly, a drive parallel to the compartment wall of the 
cavity in Fig. [lja) cannot break the mirror symmetry 
of the system in the x direction. However, it impacts 
the diffusion process through the pore, in that it presses 
the Brownian particle periodically against the horizontal 
walls, twice a period. Pressed against the top and the 
bottom walls of the cavity, symmetrically placed with 
respect to the pore, the particle is less likely to escape 
the side of the cavity it is trapped in; it is more likely to 
do so as A y (t) reverses sign, twice a period. Therefore, 
one expects N{T) to develop a denser peak structure with 



T n = ~T n , n = l,2,3,..., (9) 




FIG. 3: (Color online) Harmonic mixing in a cavity. Spectral 
density of the coordinate x(t) for two different noise inten- 
sities, D, and (a) either a longitudinal or a transverse drive 
(see legend); (b) a combination of longitudinal and trans- 
verse drives. The drive parameters are A x = 0.05, A y — 0.5, 
Cl x = 0.01, and Cl y — Cl x /\/2. The harmonic mixing reso- 
nances are indicated by down-pointing vertical arrows with 
relevant indices (m, n). 



as shown in Fig. [2jR). Analogously to the situation of 
Fig. [2|a), the distribution peaks become sharpest for an 
optimal value of Cl. An explicit numerical analysis (not 
shown) confirms that such an optimal synchronization 
occurs also under the resonance condition of Eq. ^ , but 
with the difference that now ta(D) stays for the trans- 
verse diffusion time r y = y^/3D. Although the right and 
left flows through the pore are modulated in time, we 
stress that in the presence of a transverse drive A y (t) , no 
geometric SR can occur because of the persistent mirror 
symmetry of the longitudinal motion. 
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FIG. 4: (Color online) Resonant behavior of the two low- 
est order harmonic mixing components (1,0), (-1,1). For the 
sake of comparison S^wi.q), is shown also for (A x ,A y ) = 
(0.005,0), geometric SR. All simulation parameters here are 
as in Fig. [3jb). 
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IV. HARMONIC MIXING 

The spectral characterization of geometric SR, as out- 
lined in Sec. [TTJ can provide us with more insight into the 
resonant transport mechanisms at work in a partitioned 
cavity. The analysis of Ref. [5] is equivalent to taking the 
power spectral density (PSD) ofx(t), S(ui), and evaluat- 
ing the 5-like spike, (it/2)x 2 (D)5(uj — fi), corresponding 
to the harmonic component of (x(t)) with frequency fl. 
As clearly shown in Fig. |3ja), S(u>) develops a series of 
spectral spikes at 



ui n = (2n — 1) O, 



1,2,3, 



(10) 



As well known from the SR literature pQ, the even har- 
monics of the driving frequency are absent as a conse- 
quence of the x —?• —x symmetry of the 2D Langevin 
equation ([!]). On the other hand, as anticipated in the 
final paragraph of Sec. SR cannot be induced by 
transverse drives. As expected, the PSD for A x = 
and A y ^ displayed in Fig. |3|a) exhibits no resonance 
spike. 

More interesting is the situation of Fig. ^b), where 
the Brownian particle is subjected to drives in both di- 
rections. The intrinsic nonlinearity of the longitudinal 
and transverse flows discussed in Sec. IIIII A and B makes 
a mixing between them possible, a phenomenon called 
harmonic mixing (HM). 

As discussed in Ref. [2D], key ingredients for HM in a 
ID nonlinear system are: (i) nonlinearity of the driven 
process; (ii) combination of at least two harmonic drives 
with angular frequencies tti and f^! (iii) commensura- 
tion of the driving frequencies, that is, fii/f^ = p/q, 
with p ad q relative prime numbers. Under such condi- 
tions, the system response to the external drives develops 
a hierarchy of harmonics at 



m£li + nJ7 2 m, n = 0, ±1, ±2, . . 



(11) 



Dynamical symmetries peculiar to the system can lead 
to the suppressions of subsets of the harmonics U) m , n . 

The system investigated in this section, however, is two 
dimensional, which means that HM can occur for any 
ratio of the driving frequencies. In our simulations we 
employed orthogonal harmonic drives, A x (t) and A y (t), 
with incommensurate frequencies, fl x and fl y . We know 



from Sec. [Hi] that the longitudinal flows driven by the 
transverse force alone can only resonate at the even har- 
monics of VL y , namely for u) = 2n£l y , n = 1, 2, 3, . . .. This 
is an effect of the mirror symmetry of the cavity with 
respect to the horizontal axis passing through the center 
of the pore. Moreover, as mentioned above, the mirror 
symmetry of the cavity with respect to its compartment 
wall restricts the periodic components of (x(t)} to the odd 

harmonics of A x (t), i.e., uj = (2n — l)O x , n = 1, 2,3 

In conclusion, the HM spectrum of the longitudinal 
flow through the cavity pore is expected to be 



with m = ±1, ±3, ±5, . . . (odd) and n — 0, ±1, ±2, . . ., 
with no commensuration condition on the ratio Q x /Q y . 
Our prediction is corroborated by the PSD. curves plot- 
ted in Fig. 3[b): il x /£l y is an irrational number and, still, 
all detectable spectral peaks could be identified by a pair 
of indices (m,n) according to Eq. (12). 

The practical implications of the HM of longitudinal 
and transverse drives is that, while the time modulation 
introduced by the harmonic signal A y (t), alone, cannot 
be picked up by the longitudinal current across a pore, 
adding a small longitudinal signal makes A y (t) detectable 
through the very same observable. We also stress that 
mixing spikes with n, m ^ are not necessarily small 
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FIG. 5: (Color online) Geometric SR for interacting particles. 
x(D) versus D for the double cavity shown in Fig.JTl for N 
identical particles repelling via the pair force of Eq. ( |13[ ) with 
a = 1 (a) and 6 (b). The drive parameters are A x — 0.05, 



mfl x + 2nfL, 



(12) 



0, and n x = 0.01. 



G 



with respect to the harmonics of the longitudinal sig- 
nal, n — 0. In fact, all PSD spikes, S(ui m . n ) manifest 
a SR dependence of their own on the noise intensity. 
For instance, in Fig. [4j for an appropriate D range, the 
HM harmonics (—1, 1) overshoots the fundamental com- 
ponent (1,0). 



V. PARTICLE INTERACTION 

To further explore the effects of the confining geome- 
try on the pore crossing mechanism, we consider now the 
role of particle interaction. We assume that N Brown- 
ian particles are randomly distributed between the two 
cavities of the system and repel one another through the 
standard vortex-vortex potential 16 , 



fi,j 



a 



(13) 



where i,j = 1,2, . . . ,N and i ^ j. The situation mod- 
eled here is recurrent in biological physics [31I6THU]. where 
constrained geometries often accommodate controllable 
concentrations of suspended particles. The exact form of 
fij is not relevant to the present discussion [21]. Here we 
address the dependence of the geometric SR on the pa- 
rameters N (concentration) and a (coupling) of a system 
of interacting particles. 

Let us consider the geometry SR setup of Sec. [Tljwith a 
longitudinal harmonic drive, A x (t), of fixed angular fre- 
quency Q x . In Fig. [5] we plot the curves x(D) for low 
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FIG. 6: (Color online) x(D) versus a for the same setup as in 
Fig. [5] with N and D fixed (see legend). The parameter a is 
the strength of the particle-particle interaction. x(D) shows 
a peak for relatively small values of a, and slowly decreases 
for laser values of a. 



a, in panel (a), and large a, in panel (b), and increasing 
particle concentrations. For low a, increasing N amounts 
to reducing the effective cavity volume accessible to the 
Brownian particles. Moreover, when the particle cou- 
pling is relatively weak, their interactions only exert a 
mean field effect on the pore crossing process. As a conse- 
quence, one expects that, in the regime of low couplings, 
the SR noise intensity, -D max in Eq. ([3]), weakly depends 
on N, whereas the height of the SR peak, Eq. ([2]), dimin- 
ishes with raising N. This description interprets qual- 
itatively well the results of Fig. [5ja). An exception is 
represented by the curve with N = 1, which lies under 
the curve with N = 2 and peaks at higher D. All the re- 
maining curves are centered around the same D max , their 
maxima slowly decreasing for increasing TV, with N > 2. 
As a matter of fact, when passing from N = 1 to N = 2, 
pair repulsion clearly makes pore crossing more effective 
than for a single particle. However, as argued above, 
this effect becomes negligible when adding one particle 
at larger N. 

For large a, as shown in Fig. [5jb), the pair repulsion 
can grow so strong that pair crossings become unlikely. 
Cavity switching can only be achieved by pumping more 
noise into the system, which means that the SR peaks 
must shift to higher -D max for larger concentrations. Of 
course, for very large N, the harmonic component of 
(x(t)) is suppressed independent of the value of a. 

The a dependence of the driven flow across the pore 
is further illustrated in Fig. [6] for N = 2 and N = 3. 
The curves for low noise, say, D < D max for a = 0, are 
particularly suggestive. The noise is so small that in the 
absence of interactions pore crossing happens with a time 
constant substantially larger than a half-drive period, see 
Eq. p]). When switching on the interaction, the acces- 
sible free volume per particle shrinks, which facilitates 
the escape of single particles through the pore. This ex- 
plains the raising branch of the curves x versus a. As the 
coupling a grows stronger, particle repulsion gets global, 
with the consequence that particles tend to move in de- 
formable clusters (plastic flow [55]); at exceedingly large 
a, cavity switching stops and x, decays to zero. Accord- 
ing with the above argument, the resonant dependence 
of a; on a is less appreciable for higher noise levels. 



VI. CONCLUSIONS 

We have shown how particles suspended in a parti- 
tioned cavity can diffuse across a pore or hole in the di- 
viding wall, subjected to the combined action of thermal 
fluctuations and period drives. The particle flow across 
the pore is time-modulated at the drives' frequencies, 
with amplitudes that can be optimized by controlling 
the temperature of the system. This is a geometric effect 
where the resonance condition depends on the shape of 
the cavity and on the interactions among the particles it 
contains. 

The mechanism discussed in this paper clearly does 
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not depend on the dimensionality of the cavity (experi- 
ments, e.g., on colloidal systems, can more easily be car- 
ried out in 3D geometries), but can be affected by other 
competing effects: (i) Pore structure; For a finite-size 
particle, say, a translocating molecule, the actual cross- 
ing time varies with the wall structure inside the pore 
and in the vicinity of its opening [3] ; (ii) Microfluidic ef- 
fects: The flow of the electrolytic suspension fluid across 
the pore generates inhomogeneous velocity and electrical 
fields which act on drift and the orientation of the driven 
particles [7]. System-specific effects (i) and (ii) can in 
principle be incorporated in our model by adding appro- 
priate potential terms to the Langevin equation (fTl). 
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